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In the present paper the gauge-invariant formalism is developed for perturbations of the brane- 
world model in which our universe is realized as a boundary of a higher-dimensional spacetime. 
For the background model in which the bulk spacetime is (n 4- m)-dimensional and has the spatial 
symmetry corresponding to the isometry group of a n-dimensional maximally symmetric space, 
gauge-invariant equations are derived for perturbations of the bulk spacetime. Further, for the case 
corresponding to the brane-world model in which m = 2 and the brane is a boundary invariant 
, under the spatial symmetry in the unperturbed background, relations between the gauge-invariant 

' variables describing the bulk perturbations and those for brane perturbations are derived from 

' Israel's junction condition under the assumption of Z2 symmetry. In particular, for the case in 

which the bulk spacetime is a constant-curvature spacetime, it is shown that the bulk perturbation 
{^JQ' equations reduce to a single hyperbolic master equation for a master variable, and that the physical 

^ I condition on the gauge-invariant variable describing the intrinsic stress perturbation of the brane 

, yields a boundary condition for the master equation through the junction condition. On the basis 

of this formalism, it is pointed out that it seems to be difficult to suppress brane perturbations 
corresponding to massive excitations for a brane motion giving a realistic expanding universe model. 
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I. INTRODUCTION 
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. Motivated by the M-theory AdS-CFT correspondence in string theories and the hierarchy problem in 

"nI^ the particle theory , brane-world models in which our universe is realized as a boundary of a higher-dimensional 
spacetime have been actively studied recently |]9|-^8||. In particular, for the case in which the bulk spacetime is 5- 
dimensional Anti-de Sitter spacetime and the brane is realized as a flat 4-dimensional spacetime, the gravitational 
interaction between matters in the brane is well-described by the standard one on scales much larger than the scale 
corresponding to the brane tension p^-p^ . 

Further, as an extension of the analysis to a dynamical situation, the embedding of Robertson- Walker universe 
models into 5-dimensional Anti-de Sitter and Anti-de Sitter-Schwarzschild spacetimes has been discussed by many 
people 1^^^. In such high-symmetry cases, although the evolution equation for the cosmic scale factor is modified 
from the standard one, our universe is still a dynamically closed system, and the difference in the evolution equation 
can be neglected when the energy density of the universe becomes much smaller than the brane tension. Thus the 
brane-world model gives a new world model consistent with the present day observations. However, if one goes 
^ beyond this lowest-level approximation, it is not clear whether the brane-world model is consistent with all available 
Ci ^ observations because our universe is not dynamically closed in this model |0 . 

One of the simplest ways to analyze this problem is to investigate the behavior of perturbations of the brane-world 
model. Since perturbations of the brane are inevitably associated with perturbations in the geometry of the bulk 
spacetime, such investigation will make clear whether or not the open nature of the universe dynamics is controllable. 
It will also make possible an observational test of the model in terms of the anisotropy of the cosmic microwave 
background. 

As the starting point of investigations in this line, in the present paper, we develop a gauge-invariant formalism for 
perturbations of the brane-world model. The basic approach is the same as that originally developed for 4-dimensional 
spacetime by Gerlach and Sengupta [|39|-^ and utilized by some people in the analysis of the interaction between a 
domain wall and gravitational waves in 4-dimensional spacetimes |42|-{44| . 

The formalism consists of two parts. The first one is a gauge-invariant formalism for perturbations in the geometry 
of the bulk spacetime. This problem has already been investigated by some people for the standard case in which the 
bulk spacetime is vacuum and maximally symmetric | p5[ |. In the present paper, taking account of the developing nature 
of the brane-world model, we extend the formalism to the case in which the bulk spacetime is (m -t- n)-dimensional 
and its unperturbed geometry has only the isometry corresponding to the maximally symmetric space of dimension 
n{n > 1). This symmetry is utilized to expand perturbations in terms of the harmonic functions on n-dimensional 
maximally symmetric space and define gauge-invariant variables. 
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The second part establishes relations between the gauge-invariant variables describing perturbations of the brane 
and those for the bulk perturbations. In this part we assume that m = 2 and the (n + l)-dimensional brane is 
invariant under the isometry group of the bulk in the unperturbed model. Thus the brane represents an expanding 
Robertson- Walker universe in general. 

The paper is organized as follows. In the next section we first classify perturbations into tensor, vector, and scalar 
types in terms of the tensorial behavior with respect to the maximally symmetric rt-dimensional spacetime. Then 
for each type we define the gauge-invariant variables describing perturbations of the bulk geometry and express the 
Einstein equations in terms of them. In Section III, after introducing an gauge- invariant variable describing the motion 
of the brane, we express Israel's junction condition corresponding to the Z2 symmetry in terms of it and the bulk 
variables. We will show that this gives expressions for the intrinsic perturbation variables for the brane in terms of 
the bulk variables, and a boundary condition on the latter in terms of the intrinsic stress perturbations of the brane. 
In Section IV we specialize the formalism to the standard brane-world model in which the bulk spacetime is vacuum. 
We reduce the perturbation equations to a single hyperbolic equation for a master variable Q in a, two dimensional 
spacetime and express the junction conditions in terms of the master variable. We will show that the condition that 
the anisotropic stress perturbation of the brane should vanish yields the Neumann and Dirichlet boundary condition 
on the master variable for the tensor and vector perturbations, respectively, while the boundary condition for the 
scalar perturbation is obtained from the condition on the entropy perturbation of the brane. The last condition 
becomes non-local with respect to time except for the cases in which the brane is vacuum or p — —p. Section V is 
devoted to summary and discussion. 



II. BULK PERTURBATION EQUATIONS 
A. Background spacetime 

In this section we consider perturbations of spacetime structure on (to -|- 7T,)-dimensional spacetime M which is 
locally written as a product 

A^"+" = AA™ X /C" 3 = (z*^). (1) 

Its unperturbed background geometry is given by the metric 

ds^ = gMNdz^'dz"" = gab{v)dy''dy^ + r\y)dal, (2) 

where the metric 

do'n = ^ij{x)dx^ dx-' (3) 

is that with a constant sectional curvature K on /C". We denote the covariant derivatives, the connection coefficients, 
and the curvature tensors for the three metrics ds^, gabdy'^dy'', and dcr^ as 

ds2 ^ Vm,T^l,Rmnls, (4) 
gabiy)dy-dy''^ Da,"ri{y),^Rabcdiy), (5) 
dal =^ DiXjk{x),Ri]ki{x) ^ K{-iik-iji~-ialjk)- (6) 

The expressions for V^nt and Rmnls in terms of the corresponding quantities for the metric gab{y)dy'^dy'' and da^ 
are given in Appendix 

From the symmetry structure of Gmn the energy- momentum tensor Tmn for the background bulk geometry has 
the structure 

fa^ = 0, = PS]. (7) 
Hence the Einstein equations for the bulk spacetime 

Gmn + ^gMN — k^Tmn (8) 
are reduced in the unperturbed background to 

Gab + ^gab = K'^Tab, (9) 

Gl = n{K^P - A). (10) 
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B. Gauge-transformation of perturbations 

For the infinitesimal gauge transformation represented in terms of the coordinates as Sz^^ = the metric 
perturbation Hmn = 5gMN transforms as 

5hMN = —£^igMN = —'^m£,n — ^N^M- (11) 
By decomposing the connection this yields 

5hal>=-DaS,b~Dbia. (12) 
5ha^ = -r^Da (%\ - A^a, (13) 



~5h,, = ~D,^, - D,^, - 2rD''r^^al^J. (14) 
Similarly, the gauge transformation of the perturbation of the energy- momentum tensor S{dT)MN, 

S{Sf)MN = -^^Tmn = -i'-VLTMN - Tml^n^'' ~ TjvlVmC^ (15) 

is "written as 

SiSf)ab = -eO^Tab - facDbe - ^^Dai^ (16) 

~5{5f)a^ = -T^bAe' - r''PDa{r-''ii), (17) 

~5{5f\, = ~eDa{r''PHj ~ P{D,^j + Djii). (18) 



C. Gauge-Invariant Perturbation Equations 



In general, each tensor with rank at most 2 on the maximally symmetric space IC^ is uniquely decomposed into 
components of the three types, scalar, vector, and tensor, and each component can be expanded in terms of harmonic 
functions of the same type WGl . 



1. Tensor perturbation 

First we consider the tensor perturbation, which can be expanded in terms of the harmonic tensors fij, 

(A + fc2)Ty = 0, (19) 

with the properties 

= 0, bjT{ ^ 0. (20) 

In the present paper we omit the index labeling the harmonics as well as the summation symbol with respect to the 
index, because expansion coefficients corresponding to different eigenvalues decouple on the maximally symmetric 
space. 

Here note that the eigenvalue fc^ is always non-negative under a boundary condition making the operator A self- 
adjoint in the space. In particular, fc^ = appears only for the flat space [K = 0) since the corresponding 
eigentensors satisfy -DfcTy — 0, which yields = D'^Dk^ij = nKTjk- Thus the eigentensors for = are constant 
tensors. In the framework of the expansion in the sense, such eigentensors should be discarded. Thus we assume 
fc^ > in the following unless otherwise stated. 

For the tensor perturbation the metric perturbation is expanded as 

hab^O, ha^^O, h,j = 2r^ HtT,j . (21) 

Since the infinitesimal gauge transformation ^ — (^°, ^') has no tensor component, it follows that Ht is gauge- invariant. 
Similarly, STmn is expanded as 

dfab = 0, sft = 0, sf; - ttt;., (22) 
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where tt is the gauge-invariant variable representing the tensor-type anisotropic stress perturbation. 

Inserting these expansions into the expression for 5Rij, we obtain the foUowing gauge-invariant perturbation equa- 
tion: 

n k'^ + 2K 
-DHt ~ -Dr ■ DHt + — ^ Ht = k^tt, (23) 

where □ = D°-Da is the d'Alembertian on the m-dimensional space A/"™. 



2. Vector Perturbation 

Divergence-free vector fields can be expanded in terms of the vector harmonic Vi defined by 

(A + A:2)V, = 0, 

Av' = 0. 

From this we can define the vector-type harmonic tensor as 

I 



A-v,), 



which has the properties 



A + k^ ~{n + l)K 



= 0, 



fc2 - (n - 1)K^ 
2k 



(24) 
(25) 



(26) 

(27) 
(28) 



and expands a vector-type perturbation of a 2nd-rank tensor. 

As in the case of tensor harmonics, the eigenvalue k^ is always non-negative and k^ = occurs only for K = 0, for 
which the harmonic vectors become constant vectors. Thus, for the same reason as in the tensor harmonics, we assume 
A:^ > in the following. One subtle point of the vector harmonics is that A:^ > does not imply k^ — (n + 1)K > 
for K > 0. Hence for k^ < {n + 1)K and K > 0, the vector-type tensor harmonics defined by ( |26| ) should vanish, 
which implies that V* is a Killing vector on 5". In this case it follows from ( |2^ ) that the eigenvalue should be given 
by P = {n- 1)K. 

The vector perturbation of the metric is expanded in terms of the vector harmonics as 



hab = 0, hai = rfaVi, hij ^ 2r'^ Ht^.^j , 

and the vector perturbation of the energy-momentum tensor as 

STab = 0, STt - rr-'V,, ST; = ttY). 

For the reason stated above, Ht and tt are not defined for the mode k'^ = [n — 1)K with K > Q. 
Since the infinitesimal gauge transformation ^ has only the vector component 

= 0, e. = rLY,, 

the expansion coefficients of the perturbation transform as 

~5fa = -rDa (-) , SHt - -L, = 0, Stt = 0. 
\ r J r 



(29) 
(30) 

(31) 
(32) 



Hence, except the mode k"^ = {n — 1)K for K > {), the vector perturbation is described by the three gauge-invariant 
variables r^, tt and 



Fa = fa + iDaHr. 



(33) 
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On the other hand, for the mode = {n — 1)K with K > 0, only the combination 



fb \ { fa 



F'^'=rDa['f\-rDA^] (34) 



is gauge-invariant. 

From the components 5(5" and 5G*- of the Einstein equations we obtain the following gauge-invariant perturbation 
equations except the mode k"^ = [n — 1)K with K > Q: 



fc2 - (n - 1)K . 

k 



r"+l 

„n—l rpa\ ,^2 



Fa = -2n'Ta, (35) 



£»„(r"-^F«) = -k'tt. (36) 



On the other hand, for the mode k^ = {n — 1)K with K > 0, the second equation does not appear and the first 
equation is written as 



3. Scalar Perturbation 



Prom the scalar harmonic functions 



(A + fc2)S = 0, (38) 
we can construct the scalar-type harmonic vectors Sj as 

Si = -^AS, (39) 

[A + fc2 - (n - l)K]Si = 0, (40) 
DS' = kS, (41) 

and the scalar-type harmonic tensors as 

1 ^ ^ „ 1 



-^DiDjS + -7y§, (42) 



n 



_ n-lk'^ - nK , 



0, D,E.l = ^'^ (43) 



[A -h fc^ - 2nK]Sij = 0. (44) 



In contrast to the vector and tensor harmonics, a constant function becomes the normalizable k = 0-mode for 
K > 0, for which and vanish identically. Since §i = implies §=const, no degeneracy occurs for the scalar- 
type harmonic vectors except for this constant mode, and k^ > {n — 1)K if k^ > 0. On the other hand, vanishes 
identically for fc^ = nK. For fc^ > this occurs only for K > Q. Since the spectrum of fc^ is given by fc^ = l{l-\-n—l)K 
with non-negative integer /, it corresponds to the / = 1 harmonics. For other modes fc^ > 2nK. 

A scalar perturbation of the metric is expanded in terms of the scalar harmonics as 



hab = fabS, hai = vfaSi, hij = 2r''{HL7ijS + HrSij), (45) 
and a scalar perturbation of the energy-momentum tensor as 

5 fab = TabS, 6f^ = rr^S, , 5f} = 5PS)S + tt^ • (46) 

In these expansions terms corresponding to Ht and tt for fc^ = nK > and those corresponding to fa,Hx,Ta and 
tt for /c^ = do not exist. 

For k'^{k^ — nK) ^ 0, under the infinitesimal gauge transformation 

^a=TaS, ^i=rLSi, (47) 
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these expansion coefficients transform as 



Sfab 



-DaTb - DbTa, 



5fa = -rDa 

6Xa = Ta, 

5Hl = -—L 
nr 

SHt = -L, 
r 



-T 



-Ta, 



6Tab = -T'^D.f^b - TacDbT" - TbcDaT- 

~5Ta = -{TabT''-PTa), 

r 

S{dP) = -T'^DaP, 
5tt = 0, 



where Xa is defined as 



= ^ (fa + l^aHT 



kV^ k 

Hence, in addition to tt we can construct 5 independent gauge-invariant quantities as 

F = Hl + -Ht + -D'^rXa, 
n r 

Fab = fab + DaXb + D^Xa, 

Sab = Tab + T^DaXc + T^D^Xc + X'^DcTab, 
^a = Ta-^{flXb-PXa), 
T. = 5P + X^DaP. 



(48) 
(49) 

(50) 
(51) 

(52) 

(53) 

(54) 

(55) 
(56) 

(57) 



(58) 

(59) 
(60) 

(61) 
(62) 



On the other hand, for the modes fc^(fc^ — nK) = 0, these become gauge-dependent if we define them by putting 
undefined variables to zero. 

From the components 5Gab, SGf, 5G\ and the tracelcss part of of the Einstein equations, we obtain the following 
four gauge-invariant perturbation equations for modes fc^(/e^ — nK) ^ 0: 



-UFab + DaDeF^ + DbD^F^ + n^i-D^Fab + DaF,b + DbF^a) 
J^a^cb + nb^ca — ^ ^acbd^ 



R + 2A] Fa. 



-DaDbF^ - 2n ( DaDbF + ^DarDbF + ^DbvDaF 



D^DdF^'^ + —D^rD^F^d 
r 



2n 



+ 



2n . 
r 

.2nDF - ?!^t+2lDr ■ DF + 2{n - l)^—^F 



-UF^ - -Dr ■ DF^ + ^F, 

r r 



2 

2" c 



gab = 2k Sab, 



1 



-Z^Mr-'F^) + rDa -F^ + 2(n - l)DaF 



2k £a. 



(63) 
(64) 
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n-l 



D'^rD^'Fab 



- ^^^^^^^D<^rD\ - in - 1)^ 
2 2r2 ^ ' r 



Fab 



!LzlDr-DF^-'-^ 
2r 



1 fc2 



2n 



F^ 



n(n-l) (n - l)(n - 2) fc2 - nJ-iT „ 
-{n - l)UF + '-Dr ■ DF - '-^ '- F = k^E, 



-— [2(n-2)F + F:]=«:V. 



(65) 
(66) 



For the exceptional case — nK > ( pq ) does not exist, and for the case fc^ = (64) and ( pq ) do not appear. 
The other equations still hold although each variable is gauge-dependent. 

Here, note that from the Bianchi identities not all of these equations are independent, and some combinations of 
them yield the energy-momentum conservation law for the bulk matter perturbation. For example, if we eliminate 
DbFj^ and F° in iM) using ^) and (El), we obtain 



— r^+^E"^ - -E + TT + — r°'F„b - PF^) = 0. 



(67) 



This is just the equation 5[V mT^^) ~ 0. Similarly, applying the same procedure to the divergence of (^), we obtain 
the equation mT^ ') = which is expressed as 



^Db [r"(E^-T„^F,'')]+^E„ 



n— E + T^DbF - PDaF 



'f^DbF^ - T^'^DaFbc) 



0. 



(68) 



Thus, naively speaking, only 171(771 — l)/2 components of ( |63| ) are independent under ( p4[ ) and (|66|), provided that 
the bulk energy-momentum conservation laws ( pTj) and ( |68|) are satisfied. However, it is in general difficult to extract 
such component explicitly. 



III. JUNCTION CONDITION 



In the brane- world model the bulk spacetime A4 has one or two boundaries, and we live in a boundary S. Hence the 
intrinsic geometry of E is determined by the continuity of the bulk metric gi^iN and is described by the induced metric 
f/^i,. The intrinsic metric g^^ determined in this way, however, is dependent on the location of the boundary E in the 
bulk spacetime even if the geometry of the bulk spacetime is given. Furthermore, in the spacetime with boundaries, 
the bulk geometry is not uniquely determined by an initial condition unless some appropriate boundary condition 
is imposed at E. Thus, in order for the brane-world model to be well formulated, we must give some additional 
prescription to determine the motion of branes and the boundary condition at the branes for the bulk geometry. 

In the brane-world models proposed so far, this prescription is obtained by assuming that the bulk spacetime with 
boundaries is obtained from a spacetime M with Z2 symmetry by identifying points connected by the corresponding 
Z2 transformation. The boundaries correspond to fixed points of the transformation in the original covering spacetime 
M. This implies that the hypersurface in M corresponding to a boundary E is in general a singular surface in the 
sense that the extrinsic curvatures if^^ of E on its two sides have the same absolute value but their signs are different. 
Such a singular spacetime is obtained when the surface has an intrinsic energy-momentum with finite surface density 
T 



As is shown by Israel |47|, this energy- momentum surface density is related to the difference of the extrinsic curvature 
on the two sides of the singular surface E. If we define K^^, in terms of the unit normal 71m to E as 

Kf,, = -V^n„ (69) 

and denote its value on the side in the direction of n*^ as if+^i/ and that on the other side as X-^i^, this relation is 
written as 

- = (ti^ - ^t5^^ , (70) 
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where the dimension of S is n + 1. In the brane- world model, if we choose the normal vector so that it is directed 
toward the inside of the bulk spacetime, i^+f^ = —K-^^ = if^. Hence the junction condition can be rewritten as 

n^T}^ ^ 2(X^ - KSi^). (71) 

Thus, when the intrinsic dynamics of matter in the brane is given, the motion of brane is constrained by this junction 
condition. 

In this section we express the perturbation of the above junction condition in terms of gauge-invariant variables. 
We consider only the case in which the unperturbed geometry of the brane is spatially homogeneous and isotropic. 
This implies the case to = 2 for the bulk spacetime, i.e., M. = x A^" locally, and the brane is represented by a 
manifold 

S = K X /C" 9 (T,a;^) = (x'^), (72) 
where /C" corresponds to the maximally symmetric space in the unperturbed background. 



A. Constraints 

The junction condition (^l|) together with the Hamiltonian constraint and the momentum constraint for the bulk 
spacetime gives relations between quantities intrinsic to the brane and the bulk energy- momentum density. First, 
from the momentum constraint 

V.(if;-ifj;) = -«:%^, (73) 

where V is the covariant derivative with respect to the induced metric g^n on S, and _L denotes the component along 
n, we obtain 

v,r; = -2T^x. (74) 

Thus when the bulk spacetime is vacuum, the intrinsic energy-momentum tensor is conserved. 
Secondly, from the Hamiltonian constraint 

- K^K"^ -R = 2K^f^±_ - 2A, (75) 

where R is the Ricci scalar of E, we obtain 

-R-— (tI^T'' - -T'^\ = 2k^T_i^ - 2A. (76) 
4 \ ^ n ) 

This implies that the expansion law of the brane universe is different from the one without the extra-dimension for 
which the relation 

(71 - \)R = -2k^T (77) 
holds if the cosmological constant is included in T^^ . 



B. Unperturbed brane motion 



In the unperturbed background the brane motion is described by the dependence of the coordinates on the 
proper time r of E, i.e., the set of functions ^"(t). We define the unit time- like vector by m° = y"". Here and 
from now on the over dot denotes the differentiation with respect to the proper time r. The unit normal to E in the 
unperturbed background is uniquely determined by u as 

ria = -eabu'', Ua = -iabn^ ■ (78) 

The extrinsic curvature is calculated as 

Krr - nkU^Oy, Kr, = 0, = -^S], (79) 
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and the unperturbed energy-momentum tensor of the brane is written as 



(80) 



Hence the junction condition is expressed as 

Dj_r _ 
r 

in-1) 



2n 

Dj_r 



(81) 
(82) 



The first of these equations imphes that the energy density of our universe is determined by the brane motion. If the 
equation of state of the cosmic matter is given, these equations determine the brane motion because represents 
the acceleration of the brane. Further, by differentiating the first equation by r and ehminating K^, we obtain 



n{p + p)- =2u''faA 



(83) 



This equation coincides with (7^) obtained from the momentum constraint. Here a denotes the value of r at the brane 
and represents the cosmic scale factor of the Robertson- Walker universe on the brane whose metric is written as 



(84) 



C. Perturbation of the Junction Condition 



The extrinsic curvature of the brane depends on the configuration of the brane as well as on the bulk geometry. If 
we denote the deviation of the brane configuration from the background one as 

5z^'' Z^' (r, x) - Zf^ + Z^n^\ (85) 

where Zj'^ is the component of Z^'^ parallel to the brane, the perturbation of the extrinsic curvature is in general 
expressed as 

5K^, = {£z„K)^, + V^V.^i + (i?i^i, - KlK^,)Z^ 

+na5f% + ]^Kbn''n''K^,. (86) 

The perturbation of the intrinsic metric of the brane also depends both on the perturbation of the bulk metric and 
on the brane configuration. To be explicit, these relations are expressed as 

5grr ^KbU^u^ -2Z'' + 2KIZ^, (87) 
5gr, = haru'' - D.Z^ + a^{Z,la^)\ (88) 
,5gy = + b,Zj + L»j Z, + 2a27y Z°- . (89) 

To proceed further, we must treat the tensor, the vector and the scalar perturbation separately. 



1. Tensor perturbation 

For the tensor perturbation the perturbation of the intrinsic metric of the brane is expanded in terms of the tensor 
harmonics as 

Sgrr = 0, Sgrt = 0, Sgij = 2a'^hTTij. (90) 

Since — for the tensor perturbation, hr is simply related to the bulk perturbation as hr — Ht- 

The perturbation of the energy-momentum tensor intrinsic to the brane is also expressed by a single expansion 
coefficient representing the anisotropic stress perturbation of the brane as 

(5T; = 0, 5TI = 0, 5T] = tttT}. (91) 
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On the other hand, the harmonic expansion of (pq) yields 



6K^ = 0, SK; = 0, SK'j = -D^HtT). 
Hence the junction condition d7|) reduces to the single equation 



(92) 



D±Ht = -—ttt- 



(93) 



In general, the anisotropic stress perturbation is not an independent dynamical variable and is expressed by other 
dynamical variables when the model is specified. In particular, in the linear perturbation framework, it is natural 
to assume that tt^ = for the tensor perturbation. In this case (|9^ ) gives a Neumann-type boundary condition for 
the wave equation of Ht obtained in Section II C 1 . Thus we obtain a well-posed system describing the evolution of 
perturbations. 



2. Vector perturbation 

For the vector perturbation the perturbation of the brane configuration is expressed in the harmonic expansion as 

= 0, Z^ = 0, Z, = aZY,. (94) 
On the other hand the intrinsic metric perturbation is expressed as 

^3tt — 0, 5gri — — a/?Vi, Sgtj — 2a?hT^ij. (95) 

Hence we obtain the relations 

/3 = -///- «(f) , (96) 

hT = HT--Z. (97) 
a 

If we construct the standard gauge-invariant variables for the intrinsic perturbation from these metric perturbation 
variables and the matter perturbation variables defined by 

5T; 5Tl = a{p + p){v- /3) V, , ST] = tttY] , (98) 

we obtain 



ag = -hT-f3^F,,, (99) 
V = v - [3 -jhr + F//. (100) 



Note that Z disappears in these expressions because it corresponds to an intrinsic diffeomorphism of the brane. On 
the other hand, in the present case the perturbation of the extrinsic curvature is expressed as 



SKI = 0, SK; = —e'^'Da — V,:, 5K] = —F^Y). (101) 



Inserting these equations into (|7l|), we obtain the following two equations: 



K'ip + p)V = re-''DJ^], (102) 



r 

K^nr = -2-F_i. (103) 
o 

The first of these gives the expression for the intrinsic perturbation variable in terms of the bulk perturbation variable. 
The second can be regarded as the boundary condition on the bulk perturbation equations in Section [I C 2 . It will 
be shown later that it gives a Dirichlet-type boundary condition when the bulk spacetime is vacuum. 
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3. Scalar perturbation 



For the scalar perturbation for which 

Z'' = Z'^S, Zj_ = Zj_§, Zi = aZSi, 
the harmonic expansion coefficients for the intrinsic metric perturbation defined by 

Sgrr = -2aS, Sgri = -afMi, 6gij = 2a'^{hLSjij + hrSij), 
are related to those for the bulk metric perturbation as 

a k 

h^ = Hr^ + ±Z+^Z^ + ^Z^, 
na a r 

Ht = Hx Z, 

a 

where 



(104) 
(105) 

(106) 
(107) 
(108) 
(109) 

(110) 



y = z'' - x^, Y± = z_i_- x±. 

Hence the intrinsic gauge-invariant variables constructed from these are related to the bulk gauge-invariant variables 
as 



^ = hL + lhT-Yag = F+ ^Y^, 
n k r 



^ = a - -j^iaag)- 



//// 



k;y^, 



where 



(Ill) 
(112) 

(113) 



In addition to these, we can construct gauge-invariant variables from the harmonic expansion of the intrinsic matter 
perturbation 

5T; = -6pS, 6T[ = a{p + p){v- p)Si , ST] = 6pS5^^ + tttS} , (114) 

as 



V = v- -hr, 



(115) 
(116) 

T = Sp-cldp. (117) 

Among these the last one represents the amplitiide of entropy perturbation of the matter. 

The perturbation of the extrinsic curvature are now expressed in terms of the gauge-invariant variables as 

SKI 

1 



pA = 6p- -p{v - 13), 



nil 



+KIY^ - Y^+\^^R + KlAY^ 



5KJ = k 



F^„ -[KI + 



5Kj = 


-D±F- 


a 




a \ 


— + - 







D±r 
r 

D±r 



Y'' + a 



2r 



S, 

Y^ 
a 

r 

Dj_r 



(118) 
(119) 



2\ 



Yi_ 



S5!- 



(120) 
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Hence the junction condition (|71j) yields the following three relations among the gauge-invariant variables for the bulk 
and the brane: 



DiF 



Dir 



2r 



n'^rfDaDbr 



D I r" 
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(121) 



\K, + \naDbF 



1 



D^F^ + 7;KF,n/ 



K 

Y 



n - 1 



'-F, 
2 



(122) 
(123) 
(124) 



These conditions have some features that are not shared by the vector and tensor perturbation. First, although the 
variables Z'^ and Z disappear as in the other types of perturbation, Y± = Z± — X± remains in the final expressions. 
This is because Y± defines the gauge-invariant amplitude of the perturbation of the brane motion, unlike Z'^ and Z, 
which correspond to intrinsic diffeomorphism of the brane. Secondly, from the last equation one finds that a condition 
on the anisotropic stress perturbation does not give any boundary condition on the bulk perturbation. Instead, it 
constrains the perturbation which cannot be simply attributed either to the intrinsic structure of the brane or to the 
bulk. 

Then where does the boundary condition comes from? We can find an answer to this question by closely inspecting 
the structure of the above equations. First, note that the gauge-invariants representing the perturbation of the 
intrinsic geometry of the brane are determined by the bulk variables through ( |lll| ) and (112). Meanwhile, (121) and 
(123) yield the expressions of th e gau ge-invariants A and V for the intrinsic matter in terms of the bulk variables. 
Inserting these expressions into (122), we obtain an expression for the amplitude of the entropy perturbation F in 
terms of the bulk variables. Like ttt, F is not a dynamical variable and should be expressed in terms of A, V and 
other intrinsic dynamical perturbat ion v ariables whose dynamics is determined when a model of the intrinsic matter 
is given. Hence we should regard (122) or an equation derived from it by eliminating the independent dynamical 
variable as the boundary condition on the bulk perturbation. This means that the boundary condition is dependent 
on the type of the intrinsic matter perturbation, e.g., adiabatic or isocurvature. In the next section we will show that 
this boundary condition becomes non-local with respect to the time coordinate of the brane. 



IV. MASTER VARIABLE 

As was shown in ^ the metric perturbation in the bulk spacetime for the tensor perturbation is described by the 
single gauge- invariant variable Ht, and it obeys a simple wave equation. Further the junction condition gives a simple 
boundary condition on it. In contrast, for the vector and the scalar perturbations, the bulk perturbation is described 
by multi-component variables and their equations have structures too complicated to be solved. Fortunately, in the 
case in which the unperturbed background of the bulk spacetime is vacuum ( and the 2-dimensional orbit space J^'^ is 
maximally symmetric for the scalar perturbation), we can find a single master variable for the bulk perturbation and 
reduce the perturbation equation to a single wave equation. In this section we analyze the structure of the junction 
condition in terms of that master variable. 



A. Vacuum Background 



We consider the case in which m = 2 in the notation of Section II A and Tmn — 0. Hence the bulk spacetime is 
(n -I- 2)-dimensional and has the isometry group corresponding to the n-dimensional maximally symmetric space in 
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the unperturbed background. In this case, from the generahzed BirkhofF theorem, the geometry of the background 
spacetime is given by either of the foUowing two famiUes of solutions: 

1) Pure product type {Dr = 0): 



A > : dS^{^/n/2A) x S'"(v/n(n - 1)/2A), (125) 
A < : AdS^{^/n/2\A\) x H'^{^n{n - 1)/2|A|), (126) 
A = 0:£;"+^'^ (127) 



2) Schwarzschild type(£'7' ^ 0): 

ds^ = -U{r)dt' + C7(0 ^ ^^^^^ 
2 M 

U(r)=K--^-\r^- (129) 
2A 

A = ^— -. (130) 
n[n + 1) 

The derivation of the solutions of the first family and their physical meaning were given by Nariai [^.^ . For the 
second family the following simple formulas hold: 



2 



2n(n - \)M 

■R = 2X+ , (131) 



□r ^ 2(n-l)M , , 

K - (Drf , 2M , , 

In particular, when the mass parameter M vanishes, the quantities on the left-hand side of these equations become 
constant, and the spacetime coincides with d^""*"^, AdS^^"^ and E^^^'^ for A > 0, < 0, and = 0, respectively. 

The background configuration of the brane in the Schwarzschild-type background geometry is determined by solving 



ll) and (p2h with (133). In particular, from (p3|), the same energy equation as in the no-extradimension case holds 



for the energy density of the brane. 



p= -n{p + p)-. (134) 



a 



In contrast, from (pl|) and (133), the decomposition [Dr) = —a + (_D^r) yields 



which is different from the standard expansion equation even in the case M = in the point that p is replaced by 
p^ . These equations form a closed system and determine p and a as functions of the intrinsic proper time r. When 
these functions are given, the embedding of the brane, (i(T),r(T)) is determined by rir) — air) and a solution of the 
equation 

dtV + U{a) 

In contrast to the Schwarzschild case, the background brane configuration becomes quite special for the pure product 



type background spacetime. In fact, since r=const in this case, it follows from ( pl| ) and (82) that p should vanish 
and is proportional to p. Since it is natural to assume p = for p — 0, the latter condition implies that the 
background brane motion is represented by a geodesic in the 2-dimensional constant curvature space Af. 
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B. Expression in terms of a master variable 



1. Tensor perturbation 

For the tensor perturbation the system is aheady described by a single variable. For completeness we recapitulate 
the equations for the tensor perturbation in the vacuum case. We need no further symmetry assumption on the 
unperturbed bulk geometry. 

The perturbation equation for the bulk is given by the homogeneous wave equation 

n fc^ + 2K 
-DHt - -Dr ■ DHt + 5 Ht ^ 0. (137) 

The junction condition gives the boundary condition 

D^Ht = - yT^T. (138) 



2. Vector perturbation 



For the vector perturbation on the vacuum bulk spacetime tt vanishes. Hence for > {n — 1)K, taking account 
of the fact that the orbit space Af is two-dimensional, (p^) implies that Fa is written in terms of a function Q as 



pa _ 



1 



(139) 



Hence the perturbation equation ( |35[ ) for Fa is expressed in terms of Q as 



r'^+^D, 



- {k^ - (n - l)K} Q 



= 0. 



The bulk perturbation equation is thus reduced to the single equation for the master variable given by 



(140) 



(141) 



where C is an integration constant, which can be set to zero by redefinition of fi. 

On the other hand, for the mode k"^ = {n — 1)K > 0, the gauge-invariant F^l^ has a single independent component 
and is expressed as 



F. 



(1) 



ab 



In terms of fi^^^ (5^ is expressed as 



This equation is easily solved to yield 



C 

r.n+1 



(142) 



(143) 



(144) 



where C is an integration constant. 

For k^ > {n — 1)K, the junction conditions (102) and (103) are expressed in terms of as 



a" 



1 



^'{p + p)V=^[k'-{n-l)K]n, 
1 



<J„ = 



a 



n-l 



a" 



(145) 
(146) 
(147) 
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The first equation gives a Dirichlet-type boundary condition on CI. The other two equations give expressions for the 
intrinsic gauge- invariant variables V and ag in terms of Q. Thus the initial value problem is well-posed for this system. 
The situation for the exceptional mode fc^ = {n — 1)K > is slightly different. For this mode we do not have 

the equation for ttt- However, this does not cause trouble because F^l'^ is explicitly given. The junction condition 
determines the only non-trivial gauge-invariant intrinsic to the brane, V, as 

K^a''+\p + p)V = C. (148) 

Here note that the momentum constraint (|7^) reduces to the conservation of T^^ in the present case and its 
perturbation gives 

It is easily checked that this equation is consistent with the above junction conditions. Thus the evolution of V is 
intrinsically determined and coincides with the no-extra-dimension case. In contrast, the evolution of cTg is determined 
only by solving the master equation, in contrast to the no-extra-dimension case in which Ug is related to V as [|4^ 

2K''^a'^{p + p)V = -[fc2 - {n-l)K]ag, (150) 

where k'^ denotes the gravitational constant on the brane. 



3. Scalar perturbation 



As shown in Section II C 3| , for the scal ar p erturbation on the vacuum background, ( |65| ) is automatically satisfied if 
the other three hold. Among the latter, (|6^ ) and (^6|) are written as 

(151) 
(152) 



F: = -2(n - 2)F, 



Here note that for the exceptional modes = and = nK > we do not have one or both of them. However, 
we can still assume that these equations hold by regarding missing equations as gauge conditions to fix the residual 
gauge freedom. 

As was shown by Mukohyama, in the case that the two-dimensional orbit space A/" is a constant curvature space, 
the general solutions to these equations are written in terms of a master variable f2 as 



F = r''-^F = ^(a + 2X)n, 

2n 

Fab - r"-2_p,fc = DaDtn - 



n ~ 1 , 



-□- 



A ^gab 



(153) 
(154) 



(see Appendix ^ for a simpler proof.) 

On the other hand, for the background geometry (128) with M = 0, (|63|) is reduced to the following equation: 



-UFab - -Dr ■ DFab 
r 



2A Fn. 



+ - 



D^r 



2{DaFcb + DbFca) + {n - 2) 



Dar 



DbT 



^cb - 



F,. 



^D^F+£^D,F- 
r r 



2)_2 



In terms of the master variable fJ, this equation is written as 

{DaDb + \gab)E{n) = 0, 

where 



E{n) 



un — Dr ■ on - 



+ (n - 2)A O 



(155) 



(156) 



(157) 
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As is shown in Appendix^ the general solution of (156) is written as 

E{n) = Cogoit,r) + Cm{t,r) + Car, 



(158) 



where Co, C'l and C2 are arbitrary constants. On the other hand, it is easy to see that the freedom in the definition 
of is expressed in terms of a solution to {DaDb + Xgab)^ = as 57 — > + cj. Since uj is again written as 
bj = CQgQ(t,r) + C[gi{t,r) + O^r with arbitrary constants Cq ~ Cj, the value of E{il) changes by the redefinition as 



-{P-nK){C^go + C[g,)~PC^r; K^O 
-k\C'ogo + C[gi) - {PO, - 2nAC7^)r; K = 



(159) 



From this we immediately see that Cq ~ C2 can be set to zero by an appropriate redefinition of for (k^—nK) ^ 0. 
On the other hand, only Co and Ci can be put to zero for fc^ = and K ^ Q, while only C2 can be put to zero 
for fc^ ~ nK > 0. In these cases, however, there still remains a residual gauge freedom in F and Fab- As is shown 
in Appendix ^, any solution Q to the homogeneous equation E{^) — can be set to zero by this residual gauge 
transformation, while the constants above that cannot be removed by the redefinition are just the gauge invariants 
for the exceptional modes. Thus the gauge equivalent classes of the solutions to the perturbed solutions form a 
one-dimensional space parametrized by C2 for the mode = and K > and a two-dimensional space parametrized 
by Co and Ci for the mode k^ = nK > 0. 

From now on we consider only modes with k^(k^ — nK) ^ 0. From the above argument, the master equation for 
these modes is always written as 



on 



-Dr-Dfl 



nK 



+ {n - 2)A 



In terms of the master variable, the junction conditions (121)~(124) are written as 

rD 



fc2 



-{avy - 

a a 
2^yL = kVt. 



K;n = K^- 
fcr+_c|pA 

a p + p 



— {p + p)V-^a 



Ik' 



""^ottt)-, 
nK ttt 



ak p + P 



(160) 

(161) 
(162) 
(163) 
(164) 



where 



1 



2a"-2 
D I r 



nK 



na^ 



A n 



= - 



1 



2a"-2 ^ 
-lfc2 



f7 + ( A'; (n - 1) 
nK 



+ {n-2)\]n 



{n~l)-Vl 



(165) 



(166) 



Here note that (163) is identical to the space component of the perturbation of the intrinsic conservation law of the 
energy-momentum tensor, V^T^^ = 0. Further the corresponding time component, which is written as 



(a>A)- 



k 



-{n-l)- 



a 

nK a 



V - n{p + p) $ 



fc2 



-ITT, 



(167) 



is obtained from the above junction conditions, as it should be. 
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As was discussed in Section [II C 3 , (161) and ( |162| ) are the equations determining the intrinsic gauge-invariants A 
and V. Hence the equation ( [L63D , or the equation for the intrinsic entropy perturbation T should be regarded as a 
boundary condition on the master variable. For ttt ~ 0, this expression is given by 



rD^ - 

r 



+ {2 + ncl)- 



, n 

rD^ - 
r 



+ |-n(l + w){2n - 2 + nw) 

, fc^ Dj_r ^ n 
-{n - 1)(1 + w)^^^n = 



nK 



rDi 



(168) 



where w = p/p. From this equation we immediately see that except for the special case in which p = — p, the junction 
condition yields a boundary condition that is non-local in time. 

In contrast, for the case p = —p, the junction condition yields a closed evolution equation for rDj_{Cl/r) or pA. 
To be precise, Sp and 6p becomes gauge invariant. Further, although V is ill-defined, the combination {p + p)V — 
{p + p){v — (3) — 5T^ / (a§,;) is well-defined and can have a non- vanishing val ue. I f we take these facts into account, the 
boundary condition for /? + p = is given by the equation obtained from (168) by the replacements = 0, w = — 1 
and r = 5p. 

Even in this case, the gauge invariants $ and \E' representing the intrinsic perturbations of the spatial curvature and 
the gravitational potential of the brane are determined only by solving the wave equation for under given initial 
data and a boundary condition. This is because we lack the relations that make the equations for intrinsic quantities 
closed in the no-extra-dimension case [M , 



k'^pA = (n - l)a^^{k^ - nK)^, 
(n - 2)$ + f = -K'^a^fc-^TTj,. 



(169) 
(170) 



Thus it may be difficult to find a natural initial condition for which the evolution law for the intrinsic perturbation 
becomes similar to the standard one. 



V. DISCUSSION 



In the present paper we have developed a gauge-invariant formalism for the perturbation of the brane- world model 
for which the background configuration has a spatial symmetry corresponding to a maximally symmetric space with 
a dimension n lower than the dimension n -I- m of the bulk spacetime. The formalism consisted of two parts. The 
first part gave a system of gauge-invariant equations for the perturbation of the bulk spacetime geometry. With 
applications to wider situations in mind, we derived the equations for generic values of n and m and for generic 
bulk matter. They give an extension of the formalism developed for the n = 2 and m = 2 case by Gerlach and 
Sengupta |39| . 

The second part was concerned with a situation specific to the brane- world model in which m = 2 and gave gauge- 
invariant equations for the junction condition corresponding to the Z2 symmetry along a brane with codimension one. 
As an immediate consequence, we have shown that, when the stress perturbation intrinsic to the brane is specified 
or expressed in terms of other intrinsic quantities, the junction condition yields a boundary condition at the brane(s) 
on the evolution equation for the bulk perturbation. 

In order to investigate the structure of the equations in more detail, we have introduced a master variable Q for the 
bulk perturbation and reduced the bulk perturbation equations to a single wave equation for in the case in which the 
bulk spacetime is vacuum. This reduction was already done by Mukohyama in the case in which the background 
geometry of the bulk spacetime is maximally symmetric. Since we were able to introduce the master potential for the 
scalar perturbation only in the case in which the two-dimensional orbit space has a constant curvature, the master 
equation we obtained is the same as that derived by Mukohyama. However, the master equation for the vector 
and tensor perturbations is more general and holds also in the case in which the background geometry is of the 
Schwarzschild black hole type. We have also given a proof different from that given by Mukohyama for the existence 
of the master potential for the scalar perturbation. 

We have also investigated the structure of the junction condition in terms of the master variable. In particular, 
we have shown that the boundary condition on the master potential obtained from the junction condition has a 
different structure depending on the type of perturbation: for the tensor and vector perturbations, the condition 
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that the anisotropic stress perturbation vanishes yields a Neumann-type and a Dirichlet-type boundary condition, 
respectively, while the boundary condition for the scalar perturbation is given by a condition on the intrinsic entropy 
perturbation and is non-local in time in general. 

Here, note that, although the master variable is used in an essential way in the analysis of the scalar perturbation, 
the introduction of the master variable is not the only way to make the problem tractable. For example, Fourier 
expansion of the original gauge-invariant variables in terms of time may also be used to make the equations simpler. 
If it works well, we can also treat the scalar perturbation in the Schwarzschild black hole type background. 

Although the main purpose of the present paper is to develop a formalism, we briefly discuss here a possible 
consequence of the formalism for the brane-world scenario. In the original Randall-Sundrum model, in which the 
brane is realized as a flat subspace in a 5-dimensional Anti-de Sitter spacetime, the bulk graviton modes which 
behave as massive particles inside the brane decouple from the massless mode. In our formalism this phenomenon is 
understood in the following way. 

Since n = 3, K = and M = in this case, in the units A = — 1, the gravitational wave in the bulk spacetime is 
described by Ht satisfying the wave equation 

-O^Ht = --dr irdrHT) + k'^Hr. (171) 
r 

Since the brane is static and located at r = 1, the boundary condition is given by drHx = 0. Under the Fourier 
expansion with respect to the time t, the mode Ht oc e^*'^* is a solution to the equation 

where y = l/r(l < y < oo) and /x^ = uj'^ — k^. If we require that the mode is normalizable in the generalized sense with 
respect to the natural metric drr oc dy/y^ which makes the right-hand side of the above wave equation self-adjoint, 
the spectrum of /i^ consists of two parts. One is the point spectrum /i^ = for which Ht is constant. The other is 
the continuous spectrum > for which Ht is proportional to y^Z2{^iy) where Z2 is a Bessel function of degree 2. 
Thus the general solution is written as 



Ht 



Jo 



iujt 



(173) 



The important point here is that the boundary condition is simply written as a relation between A and B. Hence 

the massless mode for which A = B = decouples from massive modes. If we apply the same argiiment to a dynamical 
case in which the brane is non-static and represents an expanding universe, the situation changes signiflcantly. In 
this case the boundary condition D^Ht = is expressed as a relation among A, B, and C. Hence all modes contain 
massive components. 

Of course, since the expansion rate of the present universe is small, one might expect that there is a mode in which 
the amplitude of the massive component is negligible. However, such a mode contains massive components with large 
amplitudes in the early phase of the universe due to rapid cosmic expansion. Hence, if the initial condition of the 
universe is imposed in the early universe as in the argmnent of quantum generation of perturbations, it is in general 
expected that the present day universe contains a non-negligible amount of massive gravitons. The situation is quite 
similar to the quantum particle creation due to cosmic expansion. Whether this problem is a crucial defect of the 
brane-world model, or it rather provides a new model of dark matter is a very interesting problem. 
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APPENDIX A: BACKGROUND QUANTITIES 
1. Connection coefficients 



2. Curvature tensors 

bed— bed, 
-. DoDbT i 
^ "J*- ^J' 

R'jki = [K-{Drfmjji-Sijjk). 



3. Ricci tensors 

n . 



Rub = ""Rab - -DaDbr, 

r 

Rni = 0, 



R\ 



Ur 
r 



Ur _ , K-{prf 



R = "'R- 2n — + n(n - 1) 



Gab = "^ab - -DaD^r - 

r 



4. Einstein tensors 

n{n -l)K- {Drf 



r 



Qab 



G) 

Gai 



l^j; {n-l){n-2)K-{Drf ^ 



0. 



APPENDIX B: PERTURBATIONS OF THE RICCI TENSORS OF THE BULK 

In general the perturbation of the Ricci tensor is expressed in terms of Kmn = Sqmn as 

25Rmn = -V^Vl/imat - ^M^Nh + ^M^hh^ + ^ Lh^ 

+RMLh% + RNhh^ — 2RMLNsh^^ , 

SR = -hMNR^"" + ^^V'^hMN - V^S/Mh. 

By decomposing the connection V into D and D we obtain 



25Rab = -nhab + DaD^hl + DbD^hl 
D^r 

+n {-DJlab + Dahcb + Dbhca) 

r 

+'^R^hcb + '^R^hca — '2™'Racbdh'"^ — '^^^ab 

1 - - Dhr D r 
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n — 2 

r 



-"lllo^Dahu + rDaDb { -h!l 



(DrY 



n + 2 



r 



- {Drf 


Ur 






r 




) hbi 









rDa ( I , 



(B4) 



2<5i?y = [2rD"'rDbh''^ + 2{n - l)D''rD\hab + 2rD"-D^rhab\ l^] 
D°'r ~ - D'^r - 

+ (n - l)——{D,haj + Djhai) + 2 — D^hkal,j 



^A/iy + 2 



(n-1) 



K ADrf Ur 



r" r 



2n 
r 



SR = DaDbh"^ + —D^rD^hab 



r 



D^rD^r ha. 



. ^^-DaD'h'^ + 2(n - 1) — D'ha^ 
1 - — . D^r 1 

1 - 

-Dh - n Dah - —Ah 



,K {Drf 



(B5) 



(B6) 



APPENDIX C: SCALAR MASTER VARIABLE 



In this appendix we show by a method different from th e pro of given in that Fab and F satisfying (151) and 
(152) are written in terms of the master variable Q, as in (153)and (154) if the 2-dimensional orbit space with the 
metric gab is a space J\f with a constant sectional curvature A. 

First, let Wab be a symmetric, traceless, and divergenceless tensor field on M. Let be a (time-like) Killing vector, 
which exists because M is maximally symmetric. If we put Wa = Wab^^ ^ from the divergenceless condition and the 
Killing equation, we obtain 



DaW'' = WabD'^S!' = 0. 



(CI) 
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In the same way, we obtain the conservation law for the combination Wab£ ^^c as 



1 



Here, from the traceless condition, this vector is related to Wa as 



Hence (C2) is written as 



e^'WaWb = 0, 



which implies that Wa is written as a gradient of a function W: 

Wa - DaW. 

( |Cl|) yields the Laplace equation 

aw = 0. 

Since the vector defined by rja = eabS,'' is orthogonal to and has the norm rjaT]'^ 

as 



(C2) 

(C3) 
(C4) 
(C5) 



(C6) 

-CaC"; the metric gab is written 

(C7) 



9ab = Jj{-£,a£.b + -na-nb): 

where U = — CaC- Utilizing this and the traceless condition, we obtain 

Wab^ WacSl = -^{Wa^b + ^acW^T^b) 
= -^{Waib + Wb^a - 9abW,e)- 

It is easil y ch ecked that the right-hand side of this equation is a symmetric, traceless, and divergenceless tensor if 
( |C5|) and (C6) are satisfied. 

In order to apply this formula to our problem, let us introduce the traceless tensor Zab as 



(C8) 



r'^-^Fab - Zab -{n- 2)r^~''Fgab. 

This tensor is not divergenceless: 

DbZl = nDa{r^-^F). 
In order to define a divergenceless tensor, let us introduce a variable as 

2nr''-^F = (□ + 2A)f7, 

and define Wab as 

Wab - Zab - (OaDbVL - ^p^Qab 



(C9) 
(CIO) 
(Cll) 

(C12) 



It is easy to check that Wab is a s ymmetric, traceless, and divergenceless tensor if A is constant, hence it is written in 
terms of a potential W as in (C8). 

Here, note that in the definition of 17 there exists a freedom of replacement ^ O + </> where (/> is a solution of the 
hyperbolic equation 



By this replacement Wab changes as 



(□ + 2A)0 = O. 



Wab W'ab = Wab " {DaDb^ + MQab)- 



(C13) 



(C14) 



Since is constrained by the hyperbolic equation, we can choose the initial condition of and dt<\> on an initial 
surface i=const so that W'„ = W[^ = 0, where t and r are the coordinates used in ( |128| ). This condition is written 
in terms of the potential W' for W'^^ as dtW — drW' — 0. For any boundary condition on W' that is linear and 
gives a well-posed initial value problem, the sol utio n satis fyin g this initial condition is W'=const, which implies that 
Wab — 0- Thus Fab and F are expressed as in (153) and (154). 
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APPENDIX D: GENERAL SOLUTION OF EQUATION (156) 



In this appendix we give the genera l sol ution to ( |l5q ) on a two-dimensional maximally symmetric space. We 
work in the coordinates {t,r) used in (|l28| ). Since the general solution for the case A = is obviously given by 
E = Co + Cit + with arbitrary constants Co ^ C2, we assume A 7^ below. 

First, note that in the {t,r) coordinates the non-vanishing Christoffel symbols are given by 



U' 



1 



U' 



^ tr ~ orr' ^ tt ~ ^ ' ^ rr ^ nrr 



2U 



2U 



From this equation the (ir)-component of (156) is written as 

^ DtDrE = U^/^dr {iJ-^/^dtE 

which yields 

E ^ h{t)U''^ + h{r). 

Inserting this expression into the (tt)-component of (|156D, we obtain 



= (A A + \9u)E = 



A 



xu 



t\2 



(Dl) 



(D2) 



(D3) 



(D4) 



where the overdot and the prime denote the differentiation with respect to t and r, respectively. Since AC/ + ([/')^/4 
\K is constant, this equation is equivalent to the following two ordinary differential equations: 



/i - XKh = c, 
^rf^ + h = 



AC/1/2' 

where c is a separation constant. The general solution of the first equation is given by 



hit) 



\c±^ + cit + Co; 

+coe^* + cie 



A' = 



On the other hand, the general solution for the equation for /2 is given by 



f2{r) 



2(-A)3/2r' 



Hence, after redefinitions of constants, the general solution including the case A = is expressed as 

E = CQgo{t,r) + Cigi{t,r) + C2r, 

where 

1; A = 0,^7^0 

go(r)= { \H'r+l; A ^ 0, X = 

t- XK ^0 



It is easy to check that this satisfies the remaining (rr)-component of (156) 

= {DrDr + Xgrr)E 



02 . U' 



dr + - 

2U U 



E. 



(D5) 
(D6) 



(D7) 

(D8) 
(D9) 

(DIO) 
(Dll) 

(D12) 
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APPENDIX E: EXCEPTIONAL MODES FOR SCALAR PERTURBATION WITH K >0 



In this appendix we show that the gauge-equivalent classes of the solutions to the perturbed Einstein equations are 
parametrized by a finite number of parameters for the exceptional modes fc^(fc^ — nK) = {K > 0) of the bulk scalar 
perturbation on a maximally symmetric background. 

First, let us consider the mode fc^ = 0. For this mode §i and vanish, and fa and Ht are undefined. Further, 
the gauge transformation is parametrized only by T^. Hence, setting the undefined variables to zero, F and Fab are 
written as = and Fat = fab, which transform under the gauge transformation as 



~5F = -:^T^ ~5Fab 



DbTa 



(El) 



For the same reason, the equations (Q) and (66), or equivalently, (151) and (151), do not exist for the mode fc^ = 0. 
However, we can recover these equations by regarding them as the gauge-fixing conditions. Then the residual gauge 
freedom is represented by Ta satisfying the following two conditions: 



0= 5 
0= ~5 

= -UTa 



F: + 2(,^ - 2)^^ 

DbF^ - 2DaF 
n - 2 



= -2i?aT^ 



(E2) 



K 



-Dr ■ DTa + -D\DaTb 



{n - 2)^ + (2n - 1)A } Ta - ^DarDbrT^' 
n-2 



+DaDbT 



-D'^rDbT", 



where 



F = r'^-^F, Fab - r"-^Fab, fa = r'^-^Ta. 



Equation (E2) implies that T" is represented by a scalar function T as 

fa ^ gab^JfoT, 



because the orbit space J\f is two dimensional. Inserting this expression into (E3), we obtain 

e'^^Db [r^DT - nrDr ■ DT + 2(n - l)KT] = 0. 
Hence, by replacing T by T-f-const, we obtain 

r^DT - nrDr ■ DT + 2(n - l)KT = 0. 



(E3) 

(E4) 

(E5) 
(E6) 
(E7) 



Since (151) and (152) hold under the above gauge conditions, any solution of the perturbed Einstein equations is 
parametrized by Q. satisfying [DaDb 4- Xgab)E{Q) = as for the generic mode. Let the set of solutions Vl to this 
equation be Sn- Then we have an onto map <&! from Sn to the space of solutions to the perturbed Einstein equations. 
The kernel of this map is spanned by the solutions of {DaDb + Xgab)^ = 0. On the other hand, F and Fab obtained 
by setting F — SF and Fab — SFab in (El) with Ta satisfying the above gauge- fixing condition is also a solution to 
the perturbed Einstein equations belonging to the trivial gauge-equivalent class. This correspondence defines a map 
$2 from the space Sg of solutions T to (E7). Then the set 5inv of gauge-equivalence classes to the perturbed Einstein 
equations is represented as Ss/^i^^2Sg- 

Here, note that iSs/ker$i is parametrized by the solution to the equation E{fl) = C^r, and hence by the initial 
data (r2,f2) on an initial surface and the constant C2. Similarly, Sg is parametrized by the initial data {T,T) for 
(E7). Therefore, by representing the condition <I'i(Sl) = $2(2^) as a relation between these initial data (and C2), we 
can determine 5inv 

Now let us undertake this program. First, by redefining — T as T, the condition $1(17) = ^2{T) is expressed as 



(DaDb + Xgab)^ = ^acDbD^T + ebcDaD^T 



2 2(n ■ 

- {Darcbc + Dbrcac) D^T + ^ 



'-e^aD-rD^Tgab 



(E8) 
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In the {t, r)-coordinates used in (128), with the help of the equations for and T, the trace and {t, r)-component of 
this equation are written as 



U 



n\' Co 2 . 



jji/2 (u-i/2^y ^ 2UT" + |c/' - 2(n - 1)^1 r + 2{n - 1)^T. 



(E9) 
(ElO) 



These equations have a solution for (T, T) when data (fi, Q) is given. 
On the other hand, the (r, r)-component is expressed as 

un" + -un' + xn^ 2u^/^ ( u-^^^t)' + -r, 

2 \ / r 



(Ell) 



and gives a constraint on C2. In fact, inserting the expression for T obtained from the trace, we obtain the condition 
C2 = 0. This implies that the set $]^^$2'5g coincides with the set of solutions to the homogeneous equation E{fl) = 0. 
Thus C2 is a gauge- invariant and parametrizes the space 5inv 

Next we examine the mode fc^ = nK . The argument is almost the same as in the above case. Now the harmonic 
scalar and the harmonic vector are non-trivial but the harmonic tensor vanishes. Hence only Ht is undefined, 
and Xa is defined as Xa = rfa/k. The gauge transformation of F and Fab are given by 



SF 



SF„ 









+ ^'^ L 







(E12) 
(E13) 



In the present case only the equation (151) is lacking. Hence we regard this as the gauge-fixing condition. Then 
the residual gauge freedom is parametrized by L satisfying the wave equation 



□L - -Dr ■ L + I nA + 2(n - 1)^ ) Z = 0, 



(E14) 



where L — L. After the redefinition — 2L/fc ^ L, the condition <I'i(ri) = ^3{L) is represented as 



(DaDb + Xgab)n = DaDbL - 



n-1 



DarDbL + DbrDaL 



n — 1 



Dr-DL+\ - {n-lf^^) L 



r 



gab- 



(E15) 



Here $3 represents the map from the space Sl of solutions L to the set of solutions to the perturbed Einstein equations. 
The trace and the (t, r)-component of this equation are written as 



un' + Xrn + — (Coffo + Cigi) = UL' 
nr 



nXr — (n — 1) — 
r 



(E16) 
(E17) 



which have a solution for {L,L) for any data (f2, Q). On the other hand, the (r, r)-component 



un" - XrQ' + Xn = UL" - {Xr + - — -U ] L' 



(E18) 



gives the constraint Cq — Ci = 0. Thus ^j,St coincides with the space of solutions to the homogeneous equation 
EiVi) = 0, and the space Sim of the gauge-equivalence classes of solutions is parametrized by the two gauge-invariant 
constants Co and Ci. 
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